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ANALYTICAL  EVALUATION  OF  A  FLUSH-MOUNTED  HYDROPHONE  ARRAY 
RESPONSE  TO  A  MODIFIED  CORCOS  TURBULENT 
WALL  PRESSURE  SPECTRUM 


INTRODUCTION 

The  purpose  of  this  work  is  to  evaluate  analytically  the 
performance  of  an  array  of  rectangular  hydrophones  mounted  on  a 
rigid  surface  in  reducing  flow  noise,  using  a  modified  Corcos 
model,  and  to  compare  the  results  calculated  for  the  modified 
Corcos  model  with  those  obtained  for  the  original  Corcos  model. 

The  analytical  model  used  in  the  present  work  is  an  array  of 

rectangular  hydrophones  flush-mounted  on  a  rigid  surface  and 

exposed  to  turbulent  flow.  The  mathematical  development  for  this 

type  of  problem  has  been  described  by  Ko  and  Nuttall.^  An 

2 

approximate  technique  has  been  developed  by  Thompson  to  perform 
the  rapid  evaluation  of  power  spectral  density  integrals.  In  the 
present  work,  the  technique  developed  in  reference  1  has  been 
extended  to  the  evaluation  of  the  array  response  to  a  modified 
Corcos  turbulent  wall  pressure  spectrum.  Using  the  original 
Corcos  model,  the  extensive  calculation  of  the  flow  induced  noise 
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received  by  a  rectangular  hydrophone  embedded  in  a  layer  of 

elastomer  backed  by  a  rigid  surface  was  made  by  Ko  and 

Schloemer.^  Further  calculations  of  turbulent  boundary  layer 

pressure  fluctuations  received  by  an  array  of  rectangular 

hydrophones  embedded  within  the  elastomer  layer  backed  by  an 

elastic  plate  with  a  finite  thickness  were  made  by  Ko  and 
4 

Schloeme  r  . 

An  analytical  development  for  determination  of  the  frequency 
spectral  density  that  gives  the  response  of  a  hydrophone  array 
mounted  on  a  rigid  surface  to  turbulent  boundary  layer  pressure 
fluctuations  is  presented.  The  double  integral  that  represents 
the  frequency  spectral  density  is  replaced  by  the  product  of  two 
single  integrals.  The  evaluation  of  the  single  integral  is  based 
on  the  contour  integral  (residue)  method.  Numerically  integrated 
results  compare  favorably  with  the  results  calculated  using  the 
exact  expression  presented  here. 

This  report  presents  mathematical  details  for  the  evaluation 

of  the  frequency  spectral  density  and  then  provides  calculations 

of  turbulent  boundary  layer  noise  reductions,  which  are  given 

relative  to  the  noise  level  calculated  for  a  point  hydrophone. 

2  -1 

The  frequency  spectral  density  is  expressed  in  dB///vPa  Hz  .  A 
number  of  parametric  studies  to  provide  hydrophone  array  design 
guidelines  for  a  quick  estimate  of  their  performance  is 
presented.  Calculated  results  were  compared  with  those  obtained 
using  the  original  Corcos  model. 


-2- 


TR  8943 


THEORETICAL  ANALYSIS 


The  response  of  an  array  of  rectangular  hydrophones 
flush-mounted  on  a  planar,  rigid  surface  to  turbulent  boundary 
layer  pressure  fluctuations  can  be  calculated  theoretically  from 
the  frequency  spectral  density  (figures  1  and  2): 


Q(  w) 


2  n 


+  00 


JI  P<W“>  S(k* 


k  )  A(k 
y  x 


k  )  dk  dk 
y  x 


Y  ' 


(1) 


where  Q(oo)  is  the  frequency  spectral  density,  P(kx,k^,w)  is  the 

turbulent  wall  pressure  spectrum,  S(kx,ky)  is  the  rectangular 

hydrophone  function,  A(k  ,k  )  is  the  rectangular  array  function, 

x  y 

co=2 nf  is  the  angular  frequency  in  rad/s,  f  is  the  frequency  in 
Hz,  and  kx  and  k^,  are  the  wavenumbers  in  the  x-  and  y-di  rect  ions , 
respectively.  The  turbulent  wall  pressure  spectrum  used  in  the 

5 

present  study,  a  modified  Corcos  model,  is  as  follows: 
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where  k  =  w/u  is  the  convective  wavenumber,  and  u  is  the 
c  c  c 

convective  flow  speed.  The  point  power  spectrum  P(o,w)  is  given 
by 


P(o,w)  =  aQ  pQ  v* 


4 


/w  / 


(3) 


where  a  ,  c.  and  c-  are  the  constants  to  be  determined,  is  the 
o  1  2  o 

fluid  density,  and  v*  is  the  friction  velocity  which  is  defined 
as 


* 

v 


<  Vpo 


1/2 

t 


(4) 


where  t  is  the  wall  shear  stress.  Using  the  definition  of  the 
o  3 

skin  friction  coefficient,  one  may  write  the  friction  velocity  as 


v*  =  (cf/2)1/2  U, 


(5) 


where  c ^  is  the  skin  friction  coefficient.  Schlichting® 
suggested  the  following  expression  for  c^.  by  fitting  the  relation 
between  the  friction  coefficient  cf  and  the  Reynolds  number  R  , 


c^  =  0.455  ( log  Rx ) 


-2.58 


(6) 


An  array  of 
The  rectangular 


rectangular  hydrophones  is  depicted  in  figure 
hydrophone  and  array  functions  are  given  by 


2. 
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S  (  k  ,  k  ) 
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and 
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(8) 


where  L  and  L  are  the  dimensions  of  the  rectangular  hydrophone, 
x  y 

dx  and  dy  are  the  central  distances  between  adjacent  hydrophones 
of  the  array,  and  N  and  N  are  the  numbers  of  hydrophones  in  the 

X  y 

x-  and  y-di rections ,  respectively. 

If  equations  (2),  (7),  and  (8)  are  substituted  into  equation 
(1),  then  equation  (1)  is  conveniently  rewritten  as 


Q(w)  =  2 n  {  PQ(w)  Q x(w)  Qy(w)  }  , 


(9) 


where 
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Qy(«) 
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For  the  convenience  of  the  ensuing  mathematical 

manipulations,  let  us  define  new  variables 
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b  =  Ly/2  » 
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(12) 


(13) 


Note  that  x  and  y  are  used  as  dummy  variables  for  the  purpose  of 
integration . 


Then,  using  the  parameters  defined  in  equation  (13),  equation 
(10)  is  rewritten  as 


P0(»)  -  2a0p02v*4  a3e/( n2«) . 


(14) 
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RESPONSE  OF  A  POINT  HYDROPHONE 

Using  the  parameters  defined  in  equation  (13),  the  response 
of  a  point  hydrophone  is  written  as 


Q(  u) 


-f  od  4-  CO 


r 

dx  r 

dy 
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OD 
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(  15) 


To  evaluate  the  first  integral  in  equation  (15),  let  x'=x-C;  then 


I 


1 


I 


[■ 


dx ' 


(16) 


The  integrated  value  of  equation  (16)  can  be  obtained  by 
evaluating  the  following  contour  integral: 


|  J  +  J  \  f(z)  dz  =  2 Jii  (Res)2=.a  ,  (17) 

^  r  -r 

where  z=x'+iy  (figure  3)  and  (Res)z=^a  is  the  residue  of  the 

2  2  2 

function  f(z)=l/(z  +a  )  at  z=ia.  By  Jordan's  lemma,  the  first 
integral  in  equation  (17)  becomes  zero  as  R->®,  and  the  residue  of 
f(z)  at  z=ia  is  obtained  as 
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Similarly,  one  obtains  the  second  integral  in  equation  (15)  as 
follows : 
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Then,  the  frequency  spectral  density  calculated  for  a  point 
hydrophone  is  written  as 
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2 n  Po(w)  I1 
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RESPONSE  OF  AN  ARRAY  OF  FINITE  HYDROPHONES 

Using  the  parameters  defined  in  equation  (13),  the  response 
of  an  array  of  hydrophones  to  turbulent  boundary  layer  pressure 
fluctuations  is  written  as 
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To  evaluate  the  first  integral  in  equation  (22),  let  x'=x-C;  then 
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a  M 


The  last  term  in  equation  (23)  can  be  expressed  in  a  series: 
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If  equation  (25)  is  substituted  into  equation  (23),  then  equation 
(23)  is  written  as 
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The  three  integrals  shown  in  the  brace  of  equation  (26)  can  be 
evaluated  by  considering  the  general  form 


f(z) 


exp ( 2i ( z+C ) ) 

~  , r ,2  ,  2 ,  2,2  ' 
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(27) 


where  z=x'+iy  (figure  4),  8o=d|m-n|  for  the  first  integral, 
8o=a+d|m-n|  for  the  second  integral,  and  8Q=a-d|m-n|  if  a>d|m-n| 
whereas  PQ=d|m-n|-a  if  a<d|m-n|  for  the  third  integral. 

Equation  (27)  can  be  evaluated  from  the  following  contour 
integrals : 
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Since  the  first  integral  in  equation  (28)  becomes  zero  as  R-»°°  by 
Jordan's  lemma,  equation  (28)  is  written  as 
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liirf  l£+p) 

J  + 


I 


f(z)  dz  =  Re{2ni(Res)  +  ni(Res)  _  . } 

Z  —  lot  Z  —  “  c, 


-R  -K-p) 

=  2  n  |exp(-28oa)  [  a(  3C2-a2  )  {  a(  2f3Qa+3  )  cos(28oC) 

+C(28Qa+l)  sin(28QC)} 


+  C( C2-3a2) { C(2poa+l)  cos ( 2  £ ) 

-  a(28oa+3)  sin(  2^2;)  ]/A  -  8q/B  1  ,  (31) 


where  A  =  4a3  {  [  a(  3  C2-a2  )  ]  2  +  [  C(  C2-3a2  )  ]  2  }  and  B=(t2  +  cx2)2. 

If  equation  (31)  is  substituted  into  equation  (26),  then  equation 
(26)  becomes 


M— 1 

I.  =  - 5 -  ^2  (M-  |m|  )  G  (  d  |  m  |  )  ,  (32) 

a  M  A  B  m=l-M 


where 


G  (  d  |  m  |  )  = 


B  exp( -2d  |  m  |  a)  |  C  j^a(  2d  |  m  |  a+3  )  cos(2d|m|C) 

+C( 2d | m | a+1 )  sin(2d|m|C)j 
+  D  j^C(2d|m|a+l)  cos(2d|m|C) 
-a(  2d  |  m  |  a+3  )  sin( 2d | m | O 1 1  -  A  d|m| 
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-  i  |b  exp[ -2 ( a+d | m | ) a]  ^  C  ^a[ 2 ( a+d | m | ) a+3 ]  cos [ 2 ( a+d | m | ) C ] 

+  C  [  2  { a+d | m | )«+l]  sin[2(a+d|m| R]) 

+  D  2 ( a+d | m | ) a+1 ]  cos [ 2 ( a+d | m | R ] 
-a[2(a+d|m| )a+3]sin[2( a+d | m |  ) C ] j  j 


-  A(  a+d  |  m  |  )  }►  —  g  (  d  |  m  |  )  , 


(33) 


where  C  =  <x(3C2-cx2)  and  D  =  CR2-3a2) 


The  expression  for  g(d|m|)  is  given  by 


g ( d | m | ) =  i  |b  exp[ -2( a-d |m| ) a]  [c (a[ 2 ( a-d | m | ) a+3 ]  cos [ 2 ( a-d | m | ) C) 

+  C [ 2 ( a-d | m | ) a+1 ]  sin[ 2 ( a-d | m |  ) U j 
+  D  (t[2(a-d|m| )a+l]  cos ( 2 ( a-d | m | ) C ] 

-  a[ 2 ( a-d | ra | ) a+3 ]  sin [ 2 ( a-d | m | R ] ) ]  +  A  (a-d|m|)|  (34) 
for  a>  d | m | ,  and 

=  j  j6  exp[ -2  ( d  |  m  | -a )  a]  j^C  ^a[  2  (  d  |  m  | -a  )  a+3  ]  cost  2  (  d  |  m  | -a  R] 

+C[ 2 ( d | m | -a ) a+1 ]  sin [ 2 ( d | m | -a ) t] ] 

+  D  [t[2(d|m|-a)a+l ]  cos[ 2 ( d | m| -a ) Cl 
-a[ 2 ( d | m | -a ) a+3 ]  sin[  2 ( d  |  m | -a  R] )  ]  +  A  (d|m|-a)j  (35) 


g(  d  |  m  |  ) 


for  a  <  d | m | . 
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Similarly,  the  second  integral  in  equation  (22)  is  obtained  as 

N-l 

J,  =  — tVt  LH  ( N—  |  n  |  )  G  (  e  |  n  |  )  ,  (36) 

Z  2b  zpN  n=l-N 

where 

G  (  e  |  n  |  )  =  2(b-e|n|)0  -  exp( -2e  |  n  [  |3) 

+  exp(-2b|3)  cosh(2e|n|3)  (37) 


for  b  >  e | n | ,  and 


G  (  e  |  n  |  )  =  exp(-2e|n|8)  [  cosh(  2b(3)-l  ]  (38) 

for  b  <  e | n | . 


It  should  be  noted  that  the  conditions  of  a  >  d|m|  (equation 
(35))  and  b  >  e|n|  (equation  (38))  represent  the  cases  where  the 
dimensions  of  each  hydrophone  exceed  the  corresponding 
separations,  which  implies  that  adjacent  rectangular  hydrophones 
are  overlapping  in  both  the  x-  and  y-di rections . 


The  frequency  spectral  density  for  a  single  hydrophone  can  be 
calculated  using  equation  (32)  with  M=1  and  equation  (36)  with 
N=1  as  follows: 
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(I2)M=1  =  t  ^/( a2AB ) ]  [  B{(3aC+CD) 

-exp(-2aa)[  { Ca( 2aa+3 ) +DC( 2aa+l ) }  cos(2£a) 

+  {CC(2aa+l  )-Da(  2aa+3 )  }  sin(2C,a)  ]}  -  Aa  ], 

and 

(J2)N=i  =  [  Jt/(  2b2p3  )  ]  [  2bf3  +  exp  ( -2b£ )  -  1  ]. 


(  39) 


(  40) 
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NUMERICAL  CALCULATIONS  AND  DISCUSSION 

The  frequency  spectral  density  for  an  array  of  rectangular 

hydrophones  mounted  on  a  rigid  surface  can  be  obtained  by 

evaluating  equation  (1).  The  integration  of  equation  (1) 

requires  the  evaluations  of  equations  (32)  and  (36).  The 

2  -1 

frequency  spectral  density  is  expressed  in  dB////Pa  Hz  .  The 
major  results  presented  in  this  work  are  the  noise  reductions 
defined  as 

NR  =  10  log10 | [ Q( w) ]/[ Q( w) ] Q |  dB ,  (41) 

where  [Q(«)]  is  the  frequency  spectral  density  calculated  for  a 
given  hydrophone  array  configuration  and  [Q(w)J  is  that 
calculated  for  a  point  hydrophone.  The  baseline  data  used  in  the 
present  study  are  given  below: 

pQ  =  1.0  g/cm^ 

c  =  150,000  cm/s 
o 

U  =20  knots 
v*  =  0.035  U. 

The  convective  flow  velocity,  uc,  is  a  function  of  frequency 

O 

as  shown  by  Bull  ,  and  an  empirically  fitted  expression  using 
Bull's  measured  data  is  obtained  as 

uc  =  U  [  0.6  +  0.4  exp( -0 . 8<*>S*/U )  ]  ,  (42) 
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•Jf 

where  6  =  0.35  cm  is  the  average  value  of  the  displacement 

thicknesses  shown  in  Bull's  measurements.  As  can  be  seen  in 
equation  (42),  the  convective  flow  velocity  asymptotically 
approaches  the  free  stream  velocity,  U,  as  the  frequency 
decreases,  and  0.6  U  as  the  frequency  increases. 

q 

A  recent  work  by  Sherman,  Ko,  and  Buehler  has  shown  that  the 
difference  between  the  peak  level  (convective  ridge)  of  the 
turbulent  wall  pressure  spectrum  and  the  level  at  the  low 
wavenumber  region  is  approximately  35  to  45  dB.  This  value  has 
been  estimated  by  adjusting  the  contamination  due  to  the  acoustic 
and  flexural  wavenumbers  at  the  low  wavenumber  region.  In  their 
study,  Martin  and  Leehey^  showed  that  their  measurements  of  low 
wavenumber  wall  pressure  spectrum  are  36  dB  below  corresponding 
convective  ridge  levels. 

Case  (1)  Modified  Corcos  Model 


Based  on  the  work  of  references  9  and  10,  a  set  of  constants 
that  makes  the  peak-to-low  wavenumber  ratio  40  dB  are  determined 
for  use  in  equation  (2)  as  follows: 

aQ  =  1,  (arbitrary  scale) 

c^  =  0.1, 

c  2  =  1.  (43) 
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The  frequency  spectral  density  has  been  evaluated  using 

equations  (32)  and  (36)  for  an  array  of  hydrophones  and  equations 

(39)  and  (40)  for  a  single  hydrophone.  Table  1  shows  the 

2  -1 

calculated  results  (dB////Pa  Hz  )  for  a  flush-mounted  point 
hydrophone  and  a  10x10  array  of  2-in. -square  hydrophones  with  no 
gaps  between  adjacent  hydrophones.  The  first  column  is  the 
frequency  in  Hz;  the  second  is  the  point  hydrophone  result;  third 
is  the  result  calculated  using  equations  (32)  and  (36);  and  the 
fourth  is  the  numerically  integrated  results  via  Simpson's  rule 
with  range  -5kc  <  (kx,k^)  <  +  5kc*  The  last  column  is  the  noise 
reduction,  which  is  the  difference  between  the  values  of  the 
second  and  third  columns. 

Figure  5  presents  the  effect  of  the  hydrophone  dimensions  on 
the  noise  reduction  as  a  function  of  frequency  for  the  flow  speed 
of  20  knots.  The  solid,  single  chain-dashed,  and  double 
chain-dashed  lines  are  the  calculated  results  for  1-,  3-  and  10- 
in. -square  hydrophones,  respectively.  The  peaks  in  the  curves 
correspond  to  locations  where  the  convective  ridge  wavenumber  is 
aligned  with  the  hydrophone  sidelobe  level  peak.  Conversely, 
nulls  in  the  noise  reduction  curves  correspond  to  the  convective 
ridge  wavenumber  aligned  with  hydrophone  beam  pattern  nulls. 
However,  it  is  shown  that  the  peaks  and  valleys  corresponding  to 
narrow  sidelobes  in  the  hydrophone  response  diminishes  as  the 
hydrophone  becomes  larger.  As  can  be  seen  in  the  figure,  more 
noise  is  attenuated  as  the  size  of  the  hydrophone  increases. 
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Figure  6  shows  the  effect  of  the  number  of  2-in. -square 
hydrophone  array  elements  on  noise  reduction  as  a  function  of 
frequency  for  the  flow  speed  of  20  knots.  These  results  are  for 
the  case  of  no  gaps  between  adjacent  hydrophones.  It  is  seen  in 
the  figure  that  as  the  number  of  array  elements  increases,  more 
noise  is  eliminated.  The  peaks  in  the  low  frequency  region 
correspond  to  locations  where  the  convective  ridge  wavenumber  is 
aligned  with  the  array  aliasing  lobe  level  peaks.  Conversely, 
nulls  in  the  noise  reduction  curves  correspond  to  the  convective 
ridge  wavenumber  aligned  with  array  beam  pattern  nulls.  However, 
it  is  shown  that  the  peaks  and  valleys  corresponding  to  narrow 
sidelobes  in  the  array  response  diminish  as  the  number  of  array 
elements  increases.  Figure  7  shows  the  results  of  a  2-in. -square 
hydrophone  array  for  the  case  of  1-in.  gaps  between  adjacent 
hydrophones . 


Case  (2)  (Original  Corcos  Model) 


8  11 

The  original  Corcos  model  is  given  as  follows:  ' 


P  (  k  ,  k  ,  co ) 

x  y 


2  4 

ao  po  v* 

co  n2 


a1  k 
1  c 


(k  -k  )2+  ( a. k  )2 
X  c  1  c 


a,k 

C 


lky+‘“2kc)2 


(44) 


where  k  =co/u  ,  and  a  ,  a.,  and  a~  are  the  constants  to  be 

C  C  OX  Zt 

determined.  Based  on  the  work  of  references  9  and  10,  a  set  of 
constants  that  makes  the  peak-to-low  wavenumber  ratio  40  dB  are 
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determined  for  use  in  equation  (43)  as  follows: 

ao  =  1,  (arbitrary  scale) 

=  0.01, 

«2  =  1.  (45) 

The  above  constants  have  been  utilized  by  Ko  and  Schloemer^ 
in  their  study  of  the  response  of  an  array  of  hydrophones 
embedded  in  an  elastomer  layer  to  a  modified  Corcos  turbulent 
wall  pressure  spectrum.  The  modified  Corcos  model  in  reference  3 
is  not  the  same  as  that  in  reference  4.  The  term  "modified"  is 
used  only  to  distinguish  the  Corcos  model  with  constants 
determined  by  Chase^.  The  constants  shown  in  reference  7  are  as 
follows : 

aQ  =  a+  ( 1+y)  -  1 

a+  =  0.766, 

r  =0.389, 

=  0.09, 

a2  =  7a^  .  (  46  ) 

It  should  be  noted  that  the  above  constants  give  a 
peak-to-low  wavenumber  ratio  of  21  dB.  These  constants  have  been 
used  by  Ko  and  Nuttall^"  in  their  study  of  the  analytical 
evaluation  of  flush-mounted  hydrophone  array  response  to  the 
Corcos  turbulent  wall  pressure  spectrum.  The  results  calculated 
using  the  constants  given  by  equation  (46)  are  presented  in  table 
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1  of  reference  1. 
reference  1  is  u 

c 


Note  that  the  convective  velocity  used  in 
=  0.6  U. 


Table  2  shows  the  results  calculated  using  the  constants 
given  by  equation  (45).  The  results  presented  in  table  2  have 
been  obtained  from  equations  (50)  and  (53)  in  reference  1  for  an 
array  of  hydrophones  and  from  equations  (28)  and  (29)  in 
reference  1  for  a  single  hydrophone,  with  the  constants  given  by 
equation  (45).  Figure  8  shows  the  effect  of  the  hydrophone 
dimensions  on  the  noise  reduction.  The  mean  flow  speed  is  20 
knots.  The  solid,  single  chain-dashed,  and  double-dashed  lines 
are  the  calculated  results  for  1-,  3-  and  10-  in. -square 
hydrophones,  respectively.  The  peaks  in  the  curves  correspond  to 
locations  where  the  convective  ridge  wavenumber  is  aligned  with 
the  hydrophone  sidelobe  level  peak.  Conversely,  nulls  in  the 
noise  reduction  curves  correspond  to  the  convective  ridge 
wavenumber  aligned  with  hydrophone  beam  pattern  nulls.  However, 
it  is  shown  that  the  peaks  and  valleys  corresponding  to  narrow 
sidelobes  in  the  hydrophone  response  diminish  as  the  hydrophone 
becomes  larger.  As  can  be  seen  in  the  figure,  more  noise  is 
attenuated  as  the  size  of  the  hydrophone  increases. 


Figure  9  shows  the  effect  of  the  number  of  2-in. -square 
hydrophone  array  elements  on  noise  reduction  as  a  function  of 
frequency  for  the  flow  speed  of  20  knots.  These  results  are  for 
the  case  of  no  gaps  between  adjacent  hydrophones.  It  is  seen  in 
the  figure  that  as  the  number  of  array  elements  increases,  more 
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noise  is  eliminated.  The  peaks  in  the  low  frequency  region 
correspond  to  locations  where  the  convective  ridge  wavenumber  is 
aligned  with  the  array  aliasing  lobe  level  peaks.  Conversely, 
nulls  in  the  noise  reduction  curves  correspond  to  the  convective 
ridge  wavenumber  aligned  with  array  beam  pattern  nulls.  However, 
it  is  shown  that  the  peaks  and  valleys  diminish  as  the  number  of 
array  elements  increases,  corresponding  to  narrow  sidelobes  in 
the  array  response.  Figure  10  shows  the  results  of  2-in. -square 
hydrophone  array  for  the  case  of  1-in.  gaps  between  adjacent 
hydrophones . 

Case  (3)  Comparison 

Figure  11  shows  a  comparison  between  the  results  calculated 
using  the  original  Corcos  model  (equation  (43)  with  constants 
given  by  equation  (45))  and  the  modified  Corcos  model  (equation 
(2)  with  constants  given  by  equation  (43))  for  a  2-in. -square 
hydrophone.  Figure  12  shows  a  comparison  between  the  results 
calculated  for  the  10x10  array  of  2-in. -square  hydrophones  with 
no  gaps  between  adjacent  hydrophones.  The  solid  and  dashed  lines 
are  the  results  calculated  for  the  original  and  modified  Corcos 
models,  respectively.  As  shown  in  the  figure,  the  results 
calculated  for  two  different  models  are  very  close  to  each  other 
for  the  low  frequency  region.  However,  as  shown  in  figures  11 
and  12,  the  results  calculated  using  the  original  Corcos  model 
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show  less  noise  reductions  compared  with  those  calculated  using 
the  modified  Corcos  model. 


In  this  study,  the  calculated  results  for  a  modified  Corcos 
model  with  constants  given  by  equation  (43)  and  the  original 
Corcos  model  with  equation  (45)  are  presented  for  comparison; 
they  are  found  to  be  in  good  agreement  as  shown  in  figures  11  and 
12  for  the  low  frequency  region.  Therefore,  the  original  Corcos 
model  adjusted  to  the  peak-to-low  wavenumber  ratio  of  40  dB  is 
recommended  in  lieu  of  the  modified  Corcos  model  for  use  in  the 
calculation  of  the  frequency  spectral  density  for  low  frequency. 
However,  the  question  of  how  to  develop  a  proper  flow  excitation 
function  still  remains  unanswered  for  the  reduction  of  turbulent 
boundary  layer  pressure  fluctuations.  Many  fundamental  questions 
regarding  extremely  complicated  turbulence  structures  are  to  be 
answered  before  an  adequate  turbulent  wall  pressure  model  is 
fully  accepted.  In  conclusion,  the  original  Corcos  model  shown 
in  equation  (44)  is  an  expression  that  shows  reasonable 
performance  in  the  low  frequency  range  when  the  empirically 
determined  constants  a^=0.01  and  0^=1  are  employed. 
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CONCLUSIONS 

An  investigation  has  been  made  for  predictions  of  the  flow 
noise  level  received  by  an  array  of  rectangular  hydrophones 
flush-mounted  on  a  rigid  surface.  The  hydrophone  array  mounted 
on  the  surface  is  exposed  to  turbulent  flow.  In  the  present 
work,  an  array  of  hydrophones  is  flush-mounted  on  a  rigid  surface 
and  is  exposed  to  turbulent  flow.  Major  results  presented  in 
this  study  are  the  turbulent  noise  reductions. 

0  The  frequency  spectral  density  that  represents  the  hydrophone 
array  response  is  obtained  by  evaluating  the  double  integral 
shown  in  equation  (1).  In  the  present  study,  the  double  integral 
is  reduced  to  the  product  of  two  single  integrals.  The  two 
single  integrals  have  been  evaluated  by  using  the  contour 
integral  method  (residue  method). 

0  It  is  found  that  more  noise  is  attenuated  as  (1)  the  number 
of  array  elements  increases,  and  (2)  the  hydrophone  dimensions 
increase.  The  array  configuration  with  no  gaps  between  adjacent 
hydrophones  provides  more  noise  attenuation  than  that  with  gaps 
between  adjacent  hydrophones  for  the  low  frequency  region. 

0  It  is  found  that  the  results  calculated  for  the  original  and 
modified  Corcos  models  are  very  close  to  each  other  when  the  wall 
pressure  spectrum  has  the  peak-to-low  wavenumber  ratio  of  40  dB . 
It  is  anticipated  that  the  original  Corcos  model  with  a^=0.01. 
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and  <*2  =  1  lieu  °f  a^  =  0.09  and  «X2=0.63)  would  be  an  acceptable 

excitation  function  for  use  in  the  flow  noise  prediction. 

0  It  should  be  remembered  that  the  peaks  and  valleys  in  noise 
reduction  curves  will  not  appear  when  a  hydrophone  or  an  array  of 
hydrophones  are  covered  by  a  layer  of  elastomer.  This  can  be 
accomplished  by  following  the  procedures  described  in  reference 
4  . 

0  This  work  provides  a  tool  for  the  quick  estimate  on  the 
response  of  an  array  of  hydrophones  to  turbulent  boundary  layer 
pressure  fluctuations.  The  analytical  expressions  shown  in  both 
the  present  work  and  reference  1  can  be  conveniently  programmed 
using  a  desk-top  computer. 

0  In  summary,  the  original  Corcos  model  shown  in  equation  (44) 
with  constants  given  by  equation  (45)  is  an  empirical  expression 
that  shows  a  reasonable  performance  when  the  frequency  spectral 
density  is  calculated.  Therefore,  the  original  Corcos  model 
adjusted  to  the  peak-to-low  ratio  of  40  dB  is  recommended  for  low 
frequency  ranges  in  lieu  of  the  modified  Corcos  model  for  use  in 
the  prediction  of  the  frequency  spectral  density. 


-25 


TR  8943 


REFERENCES 


1.  S.  H.  Ko  and  A.  H.  Nuttall,  "Analytical  Evaluation  of 
Flush-Mounted  Hydrophone  Array  Response  to  the  Corcos  Turbulent 
Wall  Pressure  Spectrum,"  Journal  of  the  Acoustical  Society  of 
America ,  Vol.  90,  No.  1  (1991) 


2.  W.  Thompson,  "The  Rapid  Evaluation  of  Two-Dimensional  Power 
Spectral  Density  Integrals  for  Large  Arrays  of  Rectangular 
Sensors  Excited  by  Turbulent  Flow,"  Applied  Research  Laboratory, 
Pennsylvania  State  University,  Technical  Memorandum  File  No. 
89-191  (29  March  1990) 

3.  S.  H.  Ko  and  H.  H.  Schloemer,  "Calculations  of  Turbulent 
Boundary  Layer  Pressure  Fluctuations  Transmitted  into  a 
Viscoelastic  Layer."  Journal  of  the  Acoustical  Society  of 
America ,  Vol.  85,  No.  4  (1989) 

4.  S.  H.  Ko  and  H.  H.  Schloemer,  "Response  of  an  Array  of 
Hydrophones  Embedded  in  an  Elastomer  Layer  to  a  Modified  Corcos 
Turbulent  Wall  Pressure  Spectrum."  Paper  presented  at  the  120th 
Meeting  of  the  Acoustical  Society  of  America,  26-30  November 
1990,  San  Diego,  CA. 

5.  R.  C.  Elswick  (Personal  Communication  1983) 


-26- 


TR  8943 


6.  H.  Schlichting,  Boundary  Layer  Theory,  McGraw-Hill,  New  York 
1968,  p.602. 

7.  D.  M.  Chase,  "The  Modelling  of  the  Wave-Vector  Frequency 
Spectrum  of  Turbulent  Boundary  Wall  Pressure,"  Journal  of  Sound 
and  Vibration,  Vol.  70  (1980) 

8.  M.  K.  Bull,  "Wall-Pressure  Fluctuations  Associated  with 
Subsonic  Turbulent  Boundary  Layer  Flow,"  Journal  of  Fluid 
Mechanics ,  Vol.  28  (1967) 

9.  C.  H.  Sherman,  S.  H.  Ko ,  and  B.  G.  Buehler,  "Measurement  of 
the  Turbulent  Boundary  Layer  of  Wave-Vector  Spectrum,"  Journal  of 
Acoustical  Society  of  America,  Vol.  88  (1990) 

10.  M.  C.  Martin  and  P.  Leehey,  "Low  Wavenumber  Wall  Pressure 
Measurements  Using  a  Rectangular  Membrane  as  a  Spatial  Filter," 
Journal  of  Sound  and  Vibration,  Vol.  52  (1977) 

11.  G.  M.  Corcos,  "The  Structure  of  the  Turbulent  Pressure  Field 
in  Boundary  Layer  Flows,"  Journal  of  Fluid  Mechanics,  Vol.  18 

(  1964  ) 


-27- 


TR  8943 


Figure  1.  Theoretical  Model 
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Configuration  of  Hydrophone  Array 
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Contour  of  Integration  (Poles  on  the  Real 
Axis)  in  the  Complex  z-Plane 


FREQUENCY  (Hz) 

Figure  5  Effect  of  the  Dimensions  of  a  Hydrophone  on  the 
Flow  Noise  Reduction  (Modified  Corcos  Model) 
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Figure  6  Effect  of  the  Number  of  Square  Hydrophone  Array 

Elements  (No  Gaps  between  Hydrophones)  on  the  Flow 
Noise  Reduction  (Modified  Corcos  Model) 
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Figure  7  Effect  of  the  Number  of  Square  Hydrophone  Array 
Elements  (Gaps  between  Hydrophones)  on  the  Flow 
Noise  Reduction  (Modified  Corcos  Model) 
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Figure  8  Effect  of  the  Dimensions  of  a  Hydrophone  on  the 
Flow  Noise  Reduction  (Original  Corcos  Model) 
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Figure  9  Effect  of  the  Number  of  Square  Hydrophone  Array 

Elements  (No  Gaps  between  Hydrophones)  on  the  Flow 
Noise  Reduction  (Original  Corcos  Model) 
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Figure  10  Effect  of  the  Number  of  Square  Hydrophone  Array 
Elements  (Gaps  between  Hydrophones)  on  the  Flow 
Noise  Reduction  (Original  Corcos  Model) 
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Figure  11  Comparison  between  the  Results  Calculated  using 
the  Original  and  Modified  Corcos  Models  for  a 
2-in. -Square  Hydrophone 
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Comparison  between  the  Results  Calculated  using 
the  Original  and  Modified  Corcos  Models  for  the 
10x10  Array  of  2-in. -Square  Hydrophones 


Figure  12 


Table  1.  Calculated  frequency  spectral  densities  for  a  10  x  10  array 
of  2-in.-square  hydrophones  with  no  gaps  between  adjacent 
hydrophones  and  for  a  point  hydrophone 
(Modified  Corcos  model  with  ci  =  0.1  and  ca=  1.0) 


Point 

Analytical 

Numerical 

Noise 

Frequency 

hydrophone 

expression 

integration 

reduction 

(Hz) 

dB//nPa2Hz'1 

dB//pPa2Hz‘1 

dB//pPa2Hz*1 

(dB) 

100 

142.5286 

102.91450 

102.91430 

-  39.34481 

200 

139.5183 

90.47374 

90.47364 

-  48.77530 

300 

137.7574 

83.08268 

83.08262 

-  54.40544 

400 

136.5080 

77.78826 

77.78821 

-  58.45048 

500 

135.5389 

73.68091 

73.68087 

-  61.58873 

600 

134.7471 

70.34100 

70.34097 

-  64.13683 

700 

134.0777 

67.54104 

67.54101 

-  66.26732 

800 

133.4977 

65.14251 

65.14249 

-  68.08593 

900 

132.9862 

63.05415 

63.05413 

-  69.66277 

1000 

132.5286 

61.21226 

61.21224 

-  71.04708 

1100 

132.1147 

59.57048 

59.57047 

-  72.27493 

1200 

131.7368 

58.09392 

58.09391 

-  73.37360 

1300 

131.3892 

56.75563 

56.75562 

-  74.36427 

1400 

131.0674 

55.53439 

55.53438 

-  75.26367 

1500 

130.7677 

54.41319 

54.41318 

-  76.08523 

1600 

130.4874 

53.37823 

53.37822 

-  76.83991 

1700 

130.2242 

52.41815 

52.41814 

-  77.53670 

1800 

129.9759 

51.52355 

51.52355 

-  78.18306 

1900 

129.7411 

50.68657 

50.68656 

-  78.78523 

2000 

129.5183 

49.90055 

49.90054 

-  79.34849 
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Table  2.  Calculated  frequency  spectral  densities  for  a  10  x  10  array 
of  2-in.-square  hydrophones  with  no  gaps  between  adjacent 
hydrophones  and  for  a  point  hydrophone 
(Original  Corcos  model  with  ai  =  0.01  and  ct2=  1.0) 


Point 

Analytical 

Numerical 

Noise 

Frequency 

hydrophone 

expression 

integration 

reduction 

(Hz) 

dB//pPa2Hz'1 

dB//(iPa2Hz'1 

dB//pPa2Hz'1 

(dB) 

100 

142.5286 

104.00510 

104.00480 

-  38.25428 

200 

139.5183 

91.42717 

91.42701 

-  47.82187 

300 

137.7574 

85.78282 

85.78272 

-  51.70530 

400 

136.5080 

80.85490 

80.85482 

-  55.38384 

500 

135.5389 

76.75074 

76.75067 

-  58.51890 

600 

134.7471 

73.91292 

73.91287 

-  60.56490 

700 

134.0777 

71.43211 

71.43207 

-  62.37624 

800 

133.4977 

69.35165 

69.35161 

-  63.87679 

900 

132.9862 

67.48494 

67.48491 

-  65.23197 

1000 

132.5286 

65.84829 

65.84826 

-  66.41105 

1100 

132.1147 

64.36728 

64.36725 

-  67.47813 

1200 

131.7368 

63.03924 

63.03921 

-  68.42829 

1300 

131.3892 

61.83305 

61.83303 

-  69.28685 

1400 

131.0674 

60.72741 

60.72739 

-  70.07065 

1500 

130.7677 

59.70817 

59.70815 

-  70.79026 

1600 

130.4874 

58.76425 

58.76423 

-  71.45389 

1700 

130.2242 

57.88625 

57.88623 

-  72.06860 

1800 

129.9759 

57.06590 

57.06589 

-  72.64071 

1900 

129.7411 

56.29619 

56.29618 

-  73.17561 

2000 

129.5183 

55.57146 

55.57144 

-  73.67758 
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